Noise and aliases in off-axis and phase-shifting holography 
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We have compared the respective efficiencies of off-axis and phase-shifting holography in terms 
of noise and aliases removal. The comparison is made by analyzing holograms of an USAF target 
backlit with laser illumination, recorded with a charge-coupled device camera. We show that it is 
essential to remove the LO beam noise, especially at low illumination levels. 
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I. INTRODUCTION 

In digital holography, holograms are recorded by a 
charge-coupled device (CCD) array detector, and the im- 
age reconstruction is performed by a computer [l|, with 
the major advantage of avoiding photographic process- 
ing. 

Off-axis holography [4] is the oldest and the simplest 
digital holography technique [3, 0, [a] . In off-axis holog- 
raphy, like in holography with photographic plates 0, 
the reference or local oscillator (LO) beam is angularly 
tilted with respect to the object observation axis. It is 
then possible to record, with a single hologram, the two 
quadratures of the object complex field. Nevertheless, 
the object field of view is reduced, since one must avoid 
the overlapping of the image with the conjugate image 
alias Q- Off-axis holography has been applied recently 
to particle Q , polarization Q , phase contrast [10| , syn- 
thetic aperture [ll| . low-coherence [12, [l3l ^^^ micro- 
scopic [iJ, [iJI imaging. 

In phase-shifting interferometry [l^ and digital holog- 
raphy [l6|, one records several images with different 
phases of the LO beam to compute the object field 
in quadrature in an on-axis (or inline) configuration. 
Recording a hologram inline requires a very accurate 
phase shift between consecutive images since the conju- 
gate image alias overlaps with the true image. Aliases are 
suppressed by making image substraction. Phase-shifting 
holography has been apphed to 3D [iT.llSl . color |19l.l20|. 
polarization [lg|, s ynt hetic aperture 21 1, low-coherence 
[22.] . surface shape [23| and microscopic i^,!^ imaging. 

Recently, we have combined the off-axis configuration 
[3, 0, M w ith our digital holography phase-shifting tech- 
nique [2l|, [2al to record off-axis phase-shifting digital 
holograms [2^ . We must notice here that our phase shift- 
ing technique minimizes phase errors 271. In the analysis 
of the data, we have used a spatial filtering technique [7] 
to remove the zero order and twin image aliases in order 
to improve the image quality. By combining all theses 



techniques, we have perform digital holography with ul- 
timate sensitivity [27] i.e. with an equivalent noise that 
reach the quantum limit of one photo electron of noise 
per reconstructed pixel during the whole measurement 
time. 

In the present paper, we propose to reaiialyze the holo- 
graphic data used in our previous paper [23 in order to 
discuss the respective merits of the different techniques 
we have combined, i.e. off-axis recording of the hologram, 
spatial filtering, 4-phase phase-shifting and low phase er- 
ror. The discussion will concern mainly aliases removal 
and sensitivity. We have chosen to consider the data of 
our previous paper [255 i^ order to simplify the discus- 
sion, since we have demonstrated in that paper that these 
data are sufficient to get optimal sensitivity. 



II. EXPERIMENTAL SET-UP 

The holographic setup is sketched in Fig[T] It con- 
sists of an interferometer in which the beams reaching the 
detector are angularly tilted (off-axis configuration) and 
frequency-shifted (dynamic phase shift). This setup has 
been described in [2^ . The main laser L is provided by a 
Sanyo DL-7140-201 diode laser (A = 780 nm, 50 mW for 
95 niA of current). It is split into an illumination beam 
(frequency lul, complex field El), and in a reference lo- 
cal oscillator (LO) beam {ujlo, Elo)- The object to be 
imaged is a back-illuminated USAF target. The optical 
frequency of the object field E is the same as the illumi- 
nation field (wi). A set of optical attenuators A (gray 
neutral filters) allows to reduce the illumination level. 
The CCD camera (PCO Pixelfiy digital camera: 12 bit, 
frame rate luqcd — 12.5 Hz, acquisition time T — 125 
ms, with 1280 x 1024 pixels of 6.7 x 6.7^m^) records the 
hologram of the object, i.e. the object {E) versus LO 
{Elo) field interference pattern. By using two acousto- 
optic modulators (AOMl and A0M2, Crystal Technol- 
ogy: uJAOMi,2 — 80 MHz), the optical frequency ojlo of 
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FIG. 1: Off-axis, dynamic phase-shifting digital holography 
setup. L: main laser; BS: Beam splitter; AOMl and AOM2; 
acousto optic modulators (Bragg cells); BE: beam expander; 
M: mirror; A: light attenuator. USAF: transmission USAF 
1951 resolution target. El, Elo and E : illumination, local 
oscillator and object field. CCD : charge-coupled device array 
detector. 



the LO beam, can be freely adjusted [27[. To make a 4- 
phase detection of the object field E, the LO frequency is 
detuned with respect to the object field optical frequency 
by: 

^L ~ l^LO — l^AOM2 — t^AOMl — ^CCd/4: (1) 

Moreover the LO beam is angularly tilted (angle 6^1°) 
with respect to the camera-to-object observation axis, so 
that the object is seen off-axis with respect to the LO 
beam propagation axis. A sequence of 12 CCD images 
/o to 111 (measurement time 0.96s) is recorded, the phase 
being dynamically shifted by tt/2 from one image to the 
next [271 , as a consequence of appropriate frequency de- 
tuning (EqlD). The interference pattern of the object 
with the LO field is carried by the frames Im recorded 
by the camera at instants t„i. 



tm = 2T:m/ujccD 



(2) 



Let us introduce the complex representations £ and 
£lo of the signal and LO field enveloppes. 



E{t) = £eJ"'* + £* 
ELo{t)^£Loe'^'+£_ 
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(3) 
(4) 



where p = —1. The grabbed CCD signal is proportional 
to the field intensity. 
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(5) 



(6) 



where c.c. is the complex conjugate of the ££1q term. 
In digital holography, the image is related to the ££1o 



term, while the zero order and twin images alias are re- 
lated to the \£lo\ and c.c. terms. The \£\ is most often 
neglected. Moreover, with the choice of the LO frequency 
lulo, the ££1q phase factor becomes 



a3'^CC0tm/4: 



jmTT/2 



J 



(7) 



III. SINGLE PHASE, OFF-AXIS 
HOLOGRAPHIC IMAGES 

In the collected data {Iq, ..., hi}, the LO beam is 
phase shifted by 7r/2 between consecutive images. To 
cancel the phase-shifting effect, we have selected images 
from the whole set of 12 recorded CCD frames for which 
the LO phase is the same : Iq, I4 and Is. A single phase, 
off-axis hologram H in the CCD plane (z — 0) is formed 
as : 



Hix, y, 0) = /o(x, y) -t- h{x, y) -t- h{x, y) 



(8) 



We have reconstructed the images by using the stan- 
dard convolution method [J, [28l that yields a calculation 
grid equal to the pixel size [23 • To calculate the con- 
volution product, we have used the Fourier method, like 
in [26l] . To avoid reconstruction aliases and to image an 
object larger than the CCD size we have enlarged our 
1280 X 1024 measurement grid by padding the data into 
a 2048 X 2048 zero matrix [30] (zero padding), as seen 
on Fig[2^, where zero is black. Further calculations are 
done onto the 2048 x 2048 grid. 

The reconstructed image is calculated by the following 
way. Eql5] yields the real space hologram H{x,y,0) in 
the CCD plane. The hologram H in the CCD reciprocal 
plane (i.e. in the z — k-space) is obtained by discrete 
Fourier transformation (FT): 



H{k^,ky,0)^FT[Hix,y,0)] 



(9) 



The k-space hologram at any distance z from the CCD 
is then: 

H{k,,ky, z) = H{k,, ky,Q) exp (jz(fc,2 + fc/)/fc) (10) 

where k — 27r/A is the optical wave vector. The 
exp [jz{kx + ky )/k) factor is the kernel function that 
describe k-space the propagation from to z. The re- 
constructed image, which is the hologram in the object 
plane (z = D), is then obtained by reverse Fourier trans- 
formation: 



Hix, y,D) = FT-' H{k^,ky,D) 



(11) 



Figl2^ shows the digital hologram H of the USAF tar- 
get recorded in the CCD plane. The hologram intensity 
\H\^ is displayed in linear gray scale. Raw CCD frames 
(1280 X 1024 pixels) are padded within a 2048 x 2048 
zero matrix displayed in black. Figl^h shows the k-space 
hologram H{kx,ky,0), whose intensity \H\'^ is displayed 
in logarithmic scale. 
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FIG. 3: Horizontal cut {y ~ 1196) of the reconstructed field 
intensity \H\'^ of Fig[2j3. Vertical axis is j_ffp in Arbitrary 
Units (A.U.). Horizontal axis is pixel index x = 0...2047. 
Vertical display is logarithmic. 



FIG. 2: Off-axis, single phase reconstruction of the tar- 
get image, (a) Hologram H{x, y, 0) that corresponds to the 
1280 X 1024 CCD image padded in a 2048 x 2048 zero matrix 
(linear gray scale display), (b) k-space field (2048 x 2048 ma- 
trix, logarithmic gray scale for the intensity |-ffp). (c) image 
calculated on the 2048 x 2048 matrix with logarithmic gray 
scale display for the intensity \H\'^. (d) Enlargement of the 
k-space true image region (dashed rectangle in FigOi). The 
400 x 400 pixels wide true image region is copied within a 
512 X 512 zero matrix. The resulting 512 x 512 \H\'^ intensity 
image is displayed in logarithmic scale. (a,d) images corre- 



spond to a total signal of : 
sequence of 3 images. 



: 4.3 X 10 photo electrons for the 



The bright zone (arrow 0) in the center of FiglJjD cor- 
responds to the FT of the LO beam (|£^loP term) that is 
the zero order image. Because the LO beam is flat in the 
CCD plane, its Fourier counterpart yields a very narrow 
bright spot located in the center of the k-space, i.e. at 

ik,,ky)C^{0,0). 

The relevant holographic signal, which corresponds to 
the S£1o interference term, is the bright zone of lower 
intensity, on the left hand side of the image (arrow -1-1). 
Because the LO beam is off-axis, the FT of the beat signal 
SE'^Q yields a signal shifted in k-space whose location can 
be precisely adjusted by tuning the angle between the LO 
and object beams (by tilting the beam splitter in front 
of the detector for example) . 

The twin image, which corresponds to the £*£lo in- 
terference term, is symmetrical to the real image with 
respect to the k-space center (arrow -1 bright zone in the 
right hand side of the image). The angular tilt (repre- 
sented by 9 on Fig[T]) between the object and the LO 
beam directions defines the separation distance between 
the true image, the twin image and the zero order image 
in k-space. 



Fig[2t shows the reconstructed image of the USAF tar- 
get whose intensity is displayed in logarithmic scale (i.e. 
\H(kx,ky, z = Up). The true image (arrow -1-1) is on 
focus for a reconstruction distance z = D = 215 mm in 
EqllOl The twin image (arrow -1) is blurred. It would 
be on focus for the reconstruction distance z = —D. The 
distance D and the off-axis angle 9 are such that the true 
and twin images are partially masked by the zero order 
alias, which is much brighter than the USAF images (±1 
orders). Nevertheless because of the logarithmic display 
the USAF true and twin images are still visible in FigO:;. 

To perform quantitative study of the Fig[2] (c) image, 
we have represented an horizontal cut along the Figl5] (c) 
dashed line. We have plotted on Fig [3] the reconstructed 
field intensity (|iif P) trace at y ~ 1196. To reduce the 
noise the curve is obtained by averaging over 11 pixels 
{y = 1191 to y = 1201). As seen on FigH:, the profile 
crosses 6 USAF black vertical bars, which are visible on 
the left hand side of the image. The profile crosses also 
the zero order image (white rectangle in the center of the 
image). On Fig[31 the zero order signal, which correspond 
to the central region of the curve {x ~ 500 to 1500), 
is much higher (about 10^ A.U.) than the USAF signal 
(about 5 X lO'^ at maximum), which is visible on the 
curve left hand side (x = to 400). The USAF bars 
(highlighted by arrows) are nevertheless visible on the 
curve. 

To select the relevant first order image, and to fully 
suppress the zero order and twin image aliases, we have 
used, as proposed by Cuche et al. Q, a k-space filtering 
(or spatial filtering) method. We have selected, in the 
k-space 2048 x 2048 matrix \H{kx,ky,z = Op, a 400 x 
400 region of interest centered on the true image bright 
zone (white dashed rectangle on Figl2]3). Note that this 
selection is made possible by the off- axis geometry that 
has translated the true image in the left hand side of the 
k-space domain. The selected area is then copied in the 
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FIG. 4: Off-axis reconstructed image of a USAF target in 
transmission with low light illumination. Images are obtained 
with k-space filtering with ~ 4.3 x 10* (a), 8.7 x 10*^ (b), 
8.7 X 10 (c) and 1.2 x 10 (d) photo electrons respectively for 



the sequence of 3 images, 
gray scale. 



|_ff| intensity displayed in linear 



center ofa512x512 zero matrix (zero padding) as shown 
on Figl2Ji- The calculation of the z = D k-space and real 
space holograms fEqfTU| are then done on this 512 x 512 
calculation grid. 

Fig [4^ shows the object plane real space hologram 
H(x, y, z) obtained by computing Eqllll which yields the 
USAF target image. Note that the translation of the se- 
lected zone in the center of the k-space domain in Figl^i 
moves the reconstructed image of the USAF target in the 
center of the image as seen in Fig|3^. The comparison 
of Figl21; with FigH^ illustrates the abihty of the spatial 
filtering method [7| to improve the quality of the recon- 
structed image in single phase off-axis digital holography. 
Nevertheless, we will see that image quality will degrade 
when the illumination level becomes low. 

The sensitivity limit of the single phase off-axis con- 
figuration was assessed by recording images of the USAF 
target at different levels of illumination. To get quanti- 
tative results, we have determined the absolute number 
of photo electrons that corresponds to the signal beam 
impinging onto the array detector. The calibration proce- 
dure is described in detail in reference [23. Figd] shows 
the reconstructed images obtained for various attenua- 
tion levels. Although the Fig[3]b, c, d images are recon- 
structed with the same holograms as in reference [25 , the 
calibration factor is slightly different. In reference [25j . 
12 holograms are used (/o ••■ In) but the holographic 
data are truncated within a 1024 x 1024 calculation ma- 
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FIG. 5: Horizontal cuts {y ~ 264) of the reconstructed field 
intensity \H\'^ corresponding to the images Fig|4] (a) to (d) 
curves correspond to (a) to (d) images. Vertical axis is \H\ 
in Arbitrary Units (A.U.). Horizontal axis is pixel index x = 
0...512. Vertical display is logarithmic. The total signal is 
~ 4.3 x 10* (a), 8.7 x lO" (b), 8.7 x 10'' (c) and 1.2 x 10"* (d) 
photo electrons respectively for the sequence of 3 images. 



trix. Here, we make use of 3 CCD frames to form the 
holograms with the whole 1280 x 1024 pixel array (cf. 
EqE]) . The number of photo electrons is thus lowered by 
a factor 1280/(4 x 1024) = 0.31 in comparison with the 
results presented in reference [25?]. 

On FigH^, with ~ 4.3 x 10^ photo electrons for all the 
pixels of the set of 3 images (i.e. for Iq + I4 + Is), one 
can see the USAF target with a good SNR (Signal to 
Noise Ratio). On Fig[4}D, with 8.7 x 10^ photo electrons, 
one can see the USAF target (arrow -1-1), but a parasitic 
signal is visible (arrow 2). When the illumination level 
goes down, the true image signal (arrow -1-1 in FigH^ and 
b) decreases, while the parasitic signal (arrow 2 in Fig|4)D, 
c and d remains unchanged. The parasitic signal (arrow 
2), which is visible on FigdJs becomes then dominant on 
Figgl:; and d with 8.7 x 10^ (c) and 1.2 x 10^ (d) photo 
electrons, while the USAF target vanishes. 

It is difficult to determine the exact nature of the par- 
asitic signal. Nevertheless, since the parasites do not 
depend on the power of signal beam, they are related to 
the LO beam alone. We can thus simply conclude that 
off-axis holography with spatial filtering is not sufficient 
to remove all LO beam parasitic contributions of our ex- 
periment, and is thus unable to reach optimal detection 
sensitivity. 

We have performed a more quantitative study of the 
off-axis images by performing cuts along the white dashed 
lines of Figd) The curves, which are obtained by by av- 
eraging over 3 pixels {y — 263 to y — 265), are displayed 
on Figini On curves (a) and (b) the signal is quite large 
(~ 4.3 X 10® and 8.7 x 10*^ photo electrons in (a) and (b)) 
and the USAF target black bars are visible (see arrows 
on Figin]). Nevertheless, because of the parasites, the 3 



left hand side bars are not visible on curve (b) . On curves 
(c) and (d) with 8.7 x 10'* and 1.2 x 10* photo electrons 
respectively, the parasitic component is dominant, and 
the USAF bars are not visible. 

The curves give the quantitative weights of the signal 
and parasitic components. Curves (b), (c) and (d) have 
roughly the same shape (except for the central region 
where the USAF bars are visible). This means that, when 
the signal is less than 8.7 x 10^ photo electrons (curve (b) 
), most of the energy lies within parasitic contributions. 



IV. PHASE-SHIFTING HOLOGRAPHIC 
IMAGES 

As mentioned above, holographic images made from 
phase- shifting measurements are reconstructed by using 
the whole set of 12 recorded CCD frames. Since the LO 
beam is phase shifted by 7r/2 between consecutive images 
(see EqlJ), the object complex hologram H' is obtained 
by summing the CCD images with the appropriate phase 
shift A(p — — m7r/2, where m = 0...11 is the image index: 



11 



H'{x,y)=Y,{-jri„,{x,y) 



(12) 



m=0 
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(13) 
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By this choice of demodulation equation (EqfT2|). the 
zero order image \£lo\ and the twin image £*£lo 
terms are both zero, since X]m=o(^-?)'" ^ ^' ^^'^ 
Em=o(-j)'"e~'^'^'^°*'"^** - 0- Moreover, the true image 
££i^o term is maximized, since (— j)™e"'"'^"^"^"*'"/'' ~ 1. 
If the 4-phase condition is not respected ujl — ^lo t^ 
ujccd/4:, the twin image term may differ from zero. This 
does not greatly affect the final result, since the zero or- 
der term, which is potentially much larger, still cancels. 

Like in single phase off-axis holography, phase-shifting 
holograms are reconstructed by the convolution method 
that involves two FTs. The 1280 x 1024 pixels wide H' 
matrix is zero padded within a square 2048 x 2048 calcu- 
lation grid and the FTs are calculated on this grid. 

^From the same data collected to assess the sensitivity 
limit of the single phase off-axis configuration, we have 
computed phase-shifting holograms of the USAF target 
for different levels of illumination, and we have recon- 
structed the images H'{x, y,z = D) of the target without 
any spatial filtering. 

Figini shows the reconstructed images whose intensity 
{\H'{x, y, -D)P) is displayed in linear scale for ~ 1.7 x 10^ 




FIG. 6: 4-phase reconstructed image of a USAF target in 
transmission with low light illumination, without k-space fil- 
tering. Images are reconstructed with 1.7 x 10^ (a), 3.5 x lO'' 
(b), 3.5 X 10^ (c) and 5 x 10* (d) photo electrons for all the 
pixels of the whole sequence of 12 images. j-ff'P intensity 
displayed in linear gray scale. 



(a), 3.5 x 10^ (b), 3.5 x 10^ (c) and 5 x lO" (d) photo 
electrons for all the pixels and all the 12 images of the 
sequence (/0.../11). These figures correspond to the same 
attenuation level that are used in the single phase off-axis 
case, but the signal, in photo electron units, is four times 
larger, since the phase-shifting holographic images are 
obtained by using 12 CCD images instead of 3 previously. 
Note that the number of pixels used in the reconstruction 
calculation is slightly larger than in reference [2^] (1280 x 
1024 pixels instead of 1024 x 1024). As a consequence, for 
the same experimental data, the total number of signal 
photo electrons must be multiplied by a factor xl.25. 

On Figin^ and b, one see the USAF target with a good 
SNR. Since the object beam is angularly tilted by 6 with 
respect to the LO beam, the USAF image is visible in the 
left hand side of the reconstructed image domain, but, 
contrarily to the Figl^t off-axis image, the true image of 
the USAF target is visible without parasitic contribution. 
Thanks to heterodyne phase-shifting [20|, the zero order 
and twin images are very low, and are thus not visible. 
On Figit, with 3.5 x 10^ photo electrons, the USAF 
target is visible, but the zero order image becomes visible 
too. On FiglBJi, with 5 x 10'' photo electrons, the SNR is 
very low. One can only guess the USAF target image in 
the left hand side of the image. 

The curves that are obtained by performing an hor- 
izontal cut {y ~ 1196) and by averaging the signal in- 
tensity |i/p over 11 pixels {y = 1191 to y = 1201), 




FIG. 7; Horizontal cuts {y ~ 1196) of the reconstructed field 
intensity \H\'^ corresponding to the 4 phase images of Fig[6] 
(a) to (d) curves correspond to (a) to (d) images. Vertical 
axis is \H\^ in Arbitrary Units (A.U.). Horizontal axis is 
pixel index x — 0...2048. Vertical display is logarithmic. The 
total signal is ~ 1.7 x 10^ (a), 3.5 x 10^ (b), 3.5 x 10^ (c) 
and 5 X 10 (d) photo electrons for all the pixels of the whole 
sequence of 12 images. 



are displayed on Fig[71 On curves (a), (b) and (c) with 
~ 1.7 X 10^, 3.5 X 10^ and 3.5 x 10^ photo electrons the 
USAF target black bars are clearly visible in the left hand 
side of the curves (see arrows on Figl?]). On curves (d) 
the signal is lower than the zero order image background 
and the bars are not visible. 

Heterodyne phase-shifting is thus very efficient for get- 
ting good images at low illumination level. Phase-shifting 
holography is clearly better in our case than off-axis 
holography with or without spatial filtering. By mak- 
ing the difference of images in Eq[l2]the |i?loP term is 
cancelled, and most of the LO beam contribution to noise 
and parasitic signal is removed by phase-shifting holog- 
raphy. This is especially important at low illumination 
level, since, in that case, the LO beam, which brings noise 
and parasites, is of much larger power than the signal. 

Since the phase-shifting holograms have been recorded 
in off-axis geometry, one can use the spatial filtering tech- 
nique to improve further the quality of the reconstructed 
images, as done in reference [2^. We have selected in the 
k-space a 400 x 400 region centered on the true image, 
copied this region in the center of a 512 x 512 calculation 
grid (zero padding) and then calculated the USAF target 
images. 

The reconstructed images are displayed on Fig|S] for 
the same levels of illumination that for Fig[6l Since the 
true image selected zone is translated in the center of the 
k-space domain, the USAF target is seen in the center of 
the reconstructed images. At high level of illumination, 
for ~ 1.7 X 10^ (a) and 3.5 x 10^ (b) photo electrons, the 
USAF images are seen with high SNR. The spatial fil- 
tering method does not seem to improve the image qual- 
ity. On Figl5](c), with 3.5 x 10^ (c) photo electrons, the 
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FIG. 8: 4-phase reconstructed image of a USAF target in 
transmission at low light illumination with k-space filtering. 
Images are reconstructed with ~ 1.7 x 10^ (a), 3.5 x lO'^ (b), 
3.5 X 10 (c) and 5 x 10 (d) photo electrons for all the pixels 



of the whole sequence of 12 images, 
in linear gray scale. 



\H'\ intensity displayed 



USAF target is still seen with a good SNR. One can 
see that the image quality is better with spatial filter- 
ing (Figl5](c)) than without (Figin](c)). Spatial filtering 
lowers the number of modes that bring noise, and cancels 
the zero order image that is seen on Fig|S] (c) (arrow 0, 
rectangular region). On Figl8](d), with 5 x lO'* (c) photo 
electrons, the USAF target is still seen with SNR ^ 1. 
As explained in reference [2^, the Fig[8](d) image is ob- 
tained with about one photo electron per resolved pixel 
of the reconstructed image (or per k-space mode), for the 
whole sequence of 12 images. 

The curves, which are obtained by averaging over 3 
pixels {y = 263 to y = 265), are displayed on FiglHl 
As seen, the USAF target black bars (see arrows) are 
easily visible on all curves. Image and cuts of FigE] and 
Figini which have been obtained with the combination of 
techniques of reference [2^ can be considered here as the 
optimal results to be reach. 



V. DISCUSSION 

Let us discuss on the respective merits of the tech- 
niques that are combined in reference [25| in order to get 
a shot-noise-limited holographic measurement. 

As seen on the Figl3] (c) reconstructed image, and on 
the Fig[3]cut, which are obtained with a large signal of 
4.3 X 10^ photo electrons, the LO beam zero order alias 




FIG. 9: Horizontal cuts {y ~ 264) of the reconstructed field 
intensity |_ffp corresponding to the images Fig|5] (a) to (d) 
curves correspond to (a) to (d) images. Vertical axis is \H\'^ 
in Arbitrary Units (A.U.). Horizontal axis is pixel index x = 
0...512. Vertical display is logarithmic. The total signal is 
~ 1.7 X 10'' (a), 3.5 X 10^ (b), 3.5 x lO'^ (c) and 5 x lO"* (d) 
photo electrons for all the pixels of the whole sequence of 12 
images. 



is much brighter than the signal itself. This is confirmed 
by FigO On the cut, the zero order alias (pixel 500 to 
1600 plateau), is about 20 times larger that the signal 
itself (peeks and valleys marked by arrows, in the pixel 
to 300 region). Moreover, the image and the zero order 
alias overlap. In single-phase regime, it is thus necessary 
to use a spatial filtering technique, as done on Fig|5] (a), 
in order to see the object without zero order alias over- 
lapping. Such spatial filtering is nevertheless far to be 
sufficient to reach the shot noise limit with our experi- 
mental data, as shown by Fig[3](c,d) and Figl5](c,d). By 
comparison, the 4-phase method without spatial filtering 
yields much better results. This is illustrated by compar- 
ing the FigH] and FigjS] images and cuts, obtained with 
1-phase holograms and spatial filtering, with the corre- 
sponding Figiniand Figl7]images and cuts, obtained with 
4-phase holograms, but without spatial filtering. 

In the case of the 1-phase images and cuts presented on 
Fig|4] and FiglSl the results are far from the shot-noise 
limit because of the parasites, visible on Fig[2] (d) (ar- 
row 2). We have done many experiments with the setup 
sketched in Fig[TJ With 1-phase, parasitic contributions 
are almost ever there, but their location depends on the 
beam splitter cube (BS) orientation. We guess that par- 
asites are related to unwanted LO beam reflections on 
the BS, whose faces are not perfectly parallel. 

To improve the quality of the 1-phase images, one can 
modify the BS orientation, such the k-space image of the 
object is moved in a "quiet zone" of the k-space, i.e. a 
zone without parasitic alias, like the upper right corner 
of the Fig[2] (b) k-space image. But this is not sufficient 
to reach the shot noise limit. To illustrate this point, we 
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FIG. 10: (a,b) k-space field intensity \Hf and \H'f recon- 
structed from 1-phase (a) and a 4-phase holograms obtained 
without signal (i.e. without illumination of the USAF char) . 
(c,d) cut along the white dashed diagonal lines for 1-phase 
(c) and 4-phase (d) holograms. Vertical axis is k-space in- 
tensity |i/p and \H'\'^ along cut. Vertical axis is pixel in- 
dex. Intensity is averaged over 11 pixels i.e. in the interval 
{i,i — 5)...(i, +5) where i = 0...2047 is the pixel index. 



have extracted, in the reference [25| experimental data, a 
sequence of 12 images without signal beam (i.e. with the 
LO beam alone). With these data, we have computed 
1-phase and 4-phase holograms (Eql8] and EafT2l respec 
tively). The k-space intensity images, \H{kj.,ky) 
\H'{kx,ky)\'^, are displayed on Fig[TO] (a) and (b) 
the same logarithmic gray scale. 

With 1-phase (a), the LO beam image is much brighter, 
and extend over a much larger area, than in the 4-phase 
case (b). To make a quantitative analysis of these im- 
ages, we have plotted profiles along the Fig[TO] (a) and 
(b) white diagonal dashed lines. Note that we have made 
diagonal cuts in order to explore " quiet zones" of the k- 
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space. The intensity signal {\H\ or |-ff'| ) along the cut 
is represented on FigUUKc) and (d). For kx = ky = 0, we 
get a peak on the (c) and (d) cuts, which corresponds to 
the flat fleld component of the LO field. In the 1-phase 
case, the peak is about ~ 4 x lO'* larger than in the 4- 
phase case. The 1-phase peak is also much broader, so 
that in most of the k-space domain, the 1-phase LO par- 
asitic signal (c) is several orders of magnitude larger than 
its 4-phase counterpart (d). It is thus much larger than 
the shot noise limit, which is equal to the 4-phase noise 
floor (within a few per cent, as verified experimentally). 
The shot noise limit is thus not reached with 1-phase de- 
tection. This means that the LO beam signal cannot be 
fully filtered-off by spatial filtering. We guess that, in real 
life experiments, it is extremely difficult to have a perfect 
flat field LO beam. Thus, complete spatial filtering the 
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FIG. 11: 1-phase reconstructed images of a USAF target in 
transmission with substraction of the average LO beam signal. 
Images are obtained with k-space filtering with ~ 4.3 x 10* 
(a), 8.7x10^ (b), 8.7x10^ (c) and 1.2x10* (d) photo electron 
respectively for the sequence of 3 images. 
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FIG. 12; Cuts of the 1-phase reconstructed images of the 
USAF target in transmission with substraction of the LO 
beam signal. Images are obtained with k-space filtering with 
~ 4.3 X 10** (a), 8.7 x 10*^ (b), 8.7 x lO" (c) and 1.2 x lO"* (d) 
photo electron respectively for the sequence of 3 images. 



statement, we have calculated, from the sequence of 12 
images, an approximation I^o of the image that should 
be obtained with the LO beam alone. Since the phase of 
the hologram interference pattern is shifted by 7r/2 from 
one CCD image to the next, the CCD signal Ilo is: 



1 " 



(14) 



It is then possible to remove the LO beam signal on each 
CCD image. The corrected 1-phase hologram H" is then 



H"{x,y,0) = Io{x,y) + h^x^y) + Is{x,y) -3lLoix,y) 

(15) 

\m=0,4,8 m=0/ 

+ J2 e+'^cc"*"/'^ ££lo + c.c. (16) 

m=0,4,8 

As seen, the \Slo\ zero order term cancels in iJ". From 
if" , we have reconstructed LO-corrected 1-phase images 
and cut, which are shown on Fig[Tl] and FiglT^l The 
attenuation level (and thus the total signal in photo elec- 
tron units) is the same than for the 1-phase images and 
cuts of FigH] and FigO Here, the images and cuts qual- 
ity is much better than without LO correction in Fig[4] 
and FigEl It is very close to the one obtained in the 
4-phase configuration in FiglH] and FiglTl With respect 
to the 4-phase case, images and cuts are slightly noisier 
here simply because the signal is smaller by a factor 4, 
since the images are reconstructed from 4 times less CCD 
frames. 

Note that these LO-corrected 1-phase results show that 
the phase accuracy in 4-phase holography is not essen- 
tial. The phase accuracy has an effect on the weight of 
the twin image alias [27] , but no significant effect on the 
weight of the zero order alias itself. To cancel the zero 
order alias, it is sufficient to build the hologram H' with 
difference of images such a way the \£\'^ term is zero. This 
is ever realized whatever the image-to-image phase shift 
is, since YlHoi-J)"" = (see EqlT3D. 



VI. CONCLUSION 



LO beam cannot be achieved. 

Performing a double filtering, in space and in time, is 
an efficient way to filter off the LO beam noise contribu- 
tions. Quite good signal-to- noise levels can be obtained 
by the time domain filtering process which consist in sub- 
tracting the LO beam signal from the 1-phase hologram 
in order to cancel the \£lo\'^ term. To illustrate this 



By comparing the reconstructed images obtained with 
holographic data measured in off-axis and phase-shifting 
configuration, we have shown that it is essential, for get- 
ting images with high SNR at low illumination levels to 
fully filter-off the LO beam. In off-axis holography, the 
LO beam yields the zero order image. It can be filtered 
off by numerical removal of its contribution in k-space 
domain (spatial filtering). This configuration also allows 



the removal of the complex conjugate image. Neverthe- 
less, this technique is not efficient enough to reach opti- 
mal sensitivity. This is particularly true when the cam- 
era exhibits image acquisition defaults. These defaults 
can be visualized during data analysis (see discussion 
about parasitic signal (arrow 2) visible on Figl2ljd)). In 
phase-shifting holography, the LO beam is filtered-off by 
making images differences. The LO beam component, 
which is the same on all the images, vanishes with im- 
ages differences. Although more efficient, phase-shifting 
is not sufficient to fully cancel the LO beam contribution 
and its noise. Single filtering (in space or time) is thus 
not enough, and double filtering (in space and time) is 
needed if optimal sensitivity is required. The key point 



of the time filtering process is to cancel the \£lo\'^ term 
in the equation yielding the hologram. We have to note 
that double filtering makes the measurement less prone 
to camera technical noise and parasites. The comparison 
of the images represented throughout this article, which 
are calculated from the same original data, is illustrative. 
Note that the quality of the phase shift in multiple-phase, 
off-axis holography is not essential, since quite as good 
results can be obtained with a single phase hologram by 
subtracting roughly the LO beam contribution. 
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